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ABSTRACT 

We have invented a novel technique to measure the radio image of a pulsar scattered by the 
interstellar plasma with 0.1 mas resolution. We extend the "secondary spectrum" analysis of 
parabolic arcs by Stinebring et al. (2001) to very long baseline interferometry and, when the 
scattering is anisotropic, we are able to map the scattered brightness astrometrically with much 
higher resolution than the diffractive limit of the interferometer. We employ this technique to 
measure an extremely anisotropic scattered image of the pulsar B0834+06 at 327 MHz. We find 
that the scattering occurs in a compact region about 420 pc from the Earth. This image has 
two components, both essentially linear and nearly parallel. The primary feature, which is about 
16 AU long and less than 0.5 AU in width, is highly inhomogeneous on spatial scales as small as 
0.05 AU. The second feature is much fainter and is displaced from the axis of the primary feature 
by about 9 AU. We find that the velocity of the scattering plasma is 16 ± 10 km s _1 approximately 
parallel to the axis of the linear feature. The origin of the observed anisotropy is unclear and 
we discuss two very different models. It could be, as has been assumed in earlier work, that the 
turbulence on spatial scales of (~ 1000 km) is homogeneous but anisotropic. However it may be 
that the turbulence on these scales is homogeneous and isotropic but the anisotropy is produced 
by highly elongated (filamentary) inhomogeneities of scale 0.05-16 AU. 



Subject headings: turbulence 
techniques: interferometric 



ISM: structure — scattering — pulsars: individual: B0834+06 



1. Introduction 

Radio pulsars provide a powerful tool for study- 
ing the ionized interstellar medium. The disper- 
sion in their pulse arrival times probes the mean 
electron density. Their very small diameters en- 
sure that they display the full range of scintilla- 
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tion and scattering phenomena, which probe the 
fine spatial structure in the electron density. Their 
pulse amplitudes exhibit a combination of diffrac- 
tive and refractive intensity scintillation on times 
from seconds to months. Compact emission from 
some bright active galactic nuclei can also show 
scintillation on times of hours to months, albeit 
smoothed by the effect of their larger angular di- 
ameters. The large body of pulsar scintillation 
data has been interpreted in terms of homogeneous 
isotropic K olmogorov turbulence in the interstel- 
lar pbasma ( Armstrong et al] 19951) . See reviews 
bv lRickettJ (|l990l k iNaravanl (|l992l ). 

However, a litany of observational evidence 
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now points to the existence of compact ionized 
structures in the interstellar medium (ISM) whose 
scattering characteristics are well beyond those 
of such homogeneous isotropic Kolmogorov tur- 
bulence. In particular the scattering is seldom 
uniformly distributed along the line of sight. It 
is often dominated by one local region somewhere 
in the line of sight, which we refer to as a "thin 
screen." Although it is unlikely to resemble a 
screen, it is thin with respect to the total line 
of sight from the source t o the observer (e.g . , 



Putney fe Stinebringl l2006t iRickett et all [2009) 



Inhomogencity in the turbulence is required on 
kiloparsec scales to explain how the level of pulsar 
scattering varies with distanc e and Galactic co- 
ordinates (jCordes et al.lll991r ). Inhomogeneity is 
also required on the parsec scale in the local ISM 
to explain the intermittent nature of the scintilla- 
tion observed in a few quasars on hour-long time 
scales (e.g., Dennett-Thorpe fc de Bruvn 2003; 
Kedziora-Chudczerll2006f )" 

In addition, evidence for AU-scale inhomogene- 
ity in the turbulence comes from extreme scatter- 
ing events (ESEs), which are observed in a few 
quasars as rare, large (10 to 50%) variations in 
flux density over s everal weeks dFiedler et alii 1987 ; 
Lazio et aij|200ll ISenkbeil et al.ll2008l ). These are 
generally seen as a decrease in flux density at- 
tributed to the passage across the lin e of sight of an 



ioniz ed cloud whi ch either scatters (|Fiedler et al 



19871 ) or refracts (|Romani et al.lll987l ) the radia- 



tion. 

In recent years there has also been increasing 
evidence that interstellar scattering (ISS) is not 
only inhomogeneous but also anisotropic. The 
most direct measure of anisotropy is through very 
long baseline interferometry (VLBI) imaging of 
scattered brightness, but it is only detectable 
on a few heavily scattered lines of sight (e.g., 



Desai fc Fevll200ll ). Anisotropy in the ISS diffrac- 



fe.g..Rickett et al. 2002 


; Dennett-Thorpe & de Bruvn 


2003 


; Coles et al. 2005; 


Bignall et al. 20061). Such 



observations give evidence for elongated fine 
structure in the ISM on scales of thousands 
of kilometers, suggesting anisotropic magneto- 
hydrodynamic turbulence co ntrolled by the mag 



netic field as discussed b y iGoldreich fc Sridhai 
(Il995l ) and lSoanglerl (|l999h . 



both homogeneity and isotropy in the ionized 
ISM. It is possible, but by no means proven, 
that the various phenomena have a common ori- 
gin in a population of AU-scale anisotropic re- 
gions of enhanced density and turbulence which 
we here generically refer to as "clouds." However, 
we note that such localized clouds must contain 
fine scale substructure that causes scattering at 
radio frequencies. Such a population of clouds 
presents a serious puzzle. Their number density 
must be many orders of magnitude greater than 
that of stars and their implied electron densities 
n c J> 10 cm -3 are much higher than expected in 
press ure equilibrium in t he warm ionized phase 



ISM ( Romani et all 1987) 



The disc overy of parabol ic arcs in the ISS of 
pulsars by IStinebring et al. ( 2001 ) has added a 



powerful new tool for probing such clouds. Many 
pulsars exhibit parabolic arcs in their secondary 
spectra (SS), which is the power spectrum (ver- 
sus delay and Doppler frequency) of the dynamic 
spectrum of intensity (versus frequency and time) . 
The arcs are sometimes narrow (in delay) which 
implies scattering by a thin layer. The distribution 
of power in the SS often reveals anisotropic scat- 
tering and in some cases there are discrete down- 
wards facing "arclets," which also imply scatter- 
ing f rom isolated anisot r opic clouds ( "c l oudle ts" ) . 
See ICordes et all (|2006h : IWalker et all (l2004h for 



interpretation of the arcs. 

The SS allows a two dimensional reconstruc- 
tion of the scattered image from observations at 
a single receiver, since each point in the SS iso- 
lates the scintillation power associated with inter- 
ference between pairs of points on the scattering 
disk. However the reconstruction can be model de- 
pendent and has an inherent two-fold amb iguity 
(jCordes et ail 120061 : iTrang fc Rickettl 120071 ) . Oc- 



casionally one can also see isolated peaks in the 



the dynamic spectrum ( 


Wolszcan & Corded 1987: 


Rickett. Lvne & Gupta 


1997). Such fringes re- 



suit from interference between the normal primary 
(on-axis) scattering disk and the off-axis discrete 
cloud. 



The examples cited above show departures from 



Recent results from lHill et al.l (|2005[ ) have com- 
pounded the difficulties in understanding the 
nature of the underlying ionized clouds. They 
observed the SS towards pulsar B0834+06 and 
found four distinct arclets scattered through 7 to 
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12 mas which they interpreted as originating from 
~0.2 AU clouds requiring n c > 100 cm -3 , similar 
to those invoked to explain ESEs towards quasars. 
By monitoring the evolution of structures in the 
secondary spectrum, they followed these clouds 
over three weeks and showed that they co-moved 
with the rest of the scattering material. 

In this paper we report VLBI observations of 
the scintillation from the same pulsar (B0834+06) 
in order to further investigate these clouds. We 
have developed a novel astrometry technique that 
makes use of SS-like quantities derived from the 
interferometer visibilities. Using these "secondary 
cross spectra" (defined in Table [T]) , we can accu- 
rately localize points on the scattering screen cor- 
responding to high signal-to-noise pixels of the SS. 
We use the results of the astrometry of many such 
points to measure the distance and velocity of the 
interstellar clouds and so define a precise model for 
the scattering. Applying this precise model to the 
astrometric image with the scattering model al- 
lows us to eliminate all the otherwise troublesome 
ambiguities and validates the model. We then use 
this precise model to recover the scattered image 
with even greater angular resolution from the sec- 
ondary spectrum itself. In our observations the 
Rayleigh resolution (i.e., synthesized beam) of the 
VLBI array is about 35 mas, the astrometric preci- 
sion is about 1 mas, and structure in the scattered 
image recovered by modeling is found on a scale 
of 100 /xas. 

2. Observations and Primary Analysis 

Pulsar B0834+06 was observed as part of a 
global VLBI project on 2005 November 12. To 
obtain the greatest astrometric precision we used 
four of the largest telescopes in the global VLBI 
network: Arecibo (AO); the Green Bank Tele- 
scope (GBT); Jodrell Bank (JB); and tied-array 
Westerbork (WB). We used a somewhat lower fre- 
quency than most secondary spectrum observa- 
tions (327 MHz) in order to obtain higher angular 
scattering and thus to better resolve the image 
with astrometry. Baseband data were recorded 
using the Mark5A disc recorders at all antennas. 
Four dual circular polarization 8 MHz bands span- 
ning the frequency range 310.5 to 342.5 MHz were 
recorded with four-level quantization, yielding a 
total data rate of 256 Mbps per antenna. In or- 



der to minimize unwanted signals that may corre- 
late between stations, the pulse calibration signals 
were disabled at all of the antennas. A total of 
5700 s of on-source data were recorded. 

An initial correlation was performed using the 
VLBA correlator in Socorro, NM with typical con- 
tinuum VLBI spectral and temporal resolutions. 
These correlator products were used for delay and 
bandpass determination. The raw data were recor- 
related at Swinbur ne University wit h the DiFX 
software correlator (jDeller et al.ll2007f ) , whose flex- 
ibility enabled the data to be processed with ex- 
tremely high spectral resolution. The output vis- 
ibilities consisting of 32768 spectral channels per 
8 MHz band (244 Hz resolution) were dumped ev- 
ery pulse period (1.25 s). The software correlator 
used incoherent dedispersion to apply a 125 ms 
wide on-pulse bin; both the on- and off-pulse vis- 
ibilities were recorded. The timing information 
from which the gate ephemeris was derived was de- 
termined through simultaneous pulsar timing ob- 
servations made at the GBT. 

Dynamic spectra were constructed from the 
autocorrelation spectra generated by the correla- 
tor for each polarization. In an analogous man- 
ner, dynamic cross spectra were formed from the 
cross correlation spectra (visibilities). In the case 
of the power spectra the off-pulse spectra were 
subtracted from the on-pulse spectra, removing 
the telescope system noise and most radio fre- 
quency interference (RFI). This was not necessary 
in forming clean dynamic cross spectra because 
the noise and RFI are not correlated between any 
two antennas. Both the dynamic spectra and the 
dynamic cross spectra were averaged into 5-pulse 
time blocks to reduce the effect of the strong pulse- 
to-pulse intensity fluctuations. The total power in 
each 8 MHz band was constrained to be constant, 
a reasonable condition considering that 8 MHz 
is much wider than the diffractive scintillation 
bandwidth, measured here as 8v ~ 3 kHz, which 
we n ote is conside rably narrower than 38 kHz 
from Cordesl ( 1986f) . This process will suppress 
any refractive intensity fluctuations, however these 
would have a time scale of many days. 

The instrumental bandpasses were corrected 
using interferometric observations of a strong 
background source. Delay, phase and amplitude 
calibration were derived from the observations of 
the target source itself using the coarsely averaged 
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visibility data from the VLBA correlator. Since 
this applies self- calibration, the astrometric posi- 
tions which we determine below are all referenced 
to the centroid of the intensity. After calibration, 
the two circular polarizations were summed to 
form the Stokes parameter / in order to maximize 
the signal to noise ratio, giving the dynamic spec- 
tra of total intensity, /(/, t), for each antenna, and 
visibility, V(f,t,b), for each baseline b. 

2.1. Secondary Spectrum Analysis and 
Arcs 

The secondary spectra, A(t, /d), are squared 
amplitudes of the two-dimensional Fourier trans- 
forms, I(t, /d), of the dynamic spectra of in- 
tensity. Here the transform variables conju- 
gate to frequency, /, and time, t, are delay, r, 
and Doppler frequency, /d, respectively. The 
quantity V(t, jb,b) represents the complex two- 
dimensional Fourier transforms of the complex 
dynamic cross spectra. As some of the quantities 
we discuss in this work have not yet been intro- 
duced in the literature, we summarize in Table[T] 
the conventions we adopt throughout for the vari- 
ous products derived from the correlator auto and 
cross correlations. 

The dynamic spectra have a spectral resolu- 
tion of 244 Hz over 8 MHz and temporal reso- 
lution of 6.25 s (about 5 pulsar rotation periods) 
over the observation duration of 6500 s, which are 
well-suited to resolving the diffractive scintillation, 
whose time scale is ~ 1 min and frequency scale 
is ~ 3 kHz. Thus the secondary spectra have a 
resolution of 125 ns in delay out to a maximum 
of 2.05 ms. They have resolution in Doppler fre- 
quency of 0.15 mHz over a width of ± 80 mHz. 
Custom software was used for this and all subse- 
quent data processing. 

For the first time, dynamic spectra of right 
minus left circular polarization and the corre- 
sponding secondary spectra were computed to 
test for differential Faraday rotation in the ISM 
(jMacquart k Melrose] 120001 ). No detectable sig- 
nal, significant at the 0.1% level, was found in any 
of these differenced spectra. In a simple model this 
implies rotation measure differences across the im- 
age of less than ~ 1.2 x 10~ 3 radm" 2 on AU scales 
(jMacquart fc Melrosdlioooh . 

The amplitude of the quantity V(t, /b,b) 



V(—T,—fr>,b) is plotted in Figure [TJ This quan- 
tity is a generalization of the secondary spectrum 
to the interferometric case and will be discussed 
later. The plot extends to a delay of 2 ms, which 
is considerably greater than has been published 
before, and is plotted with higher Doppler fre- 
quency and delay resolution than is usual. It thus 
reveals some interesting features. It is composed 
of a dense forest of fine arclets in which the apex 
of each arclet lies on or near the upwards facing 
main arc and all arcs are parabolic with th e same 
curvature, as shown by Iffill et all (|2005h . The 
disjoint group of arclets near 1 ms in delay with 
negative Doppler frequency is particularly strik- 
ing and will be referred to as the 1 ms feature. 
Their apexes lie inside the main arc and, as we 
find below, they form a separate and distinct part 
of the scattered image. Some faint arclets on the 
positive Doppler frequency side can be identified 
with delays as high as 2 ms, although they are not 
readily visible in this figure. 

3. Theory of Secondary Cross Spectra 

Here we briefly outli ne the theory of the arcs 
and refer the reade r to Walker et al. (2_004) and 
Cordes et al. (2006) for more complete analyses. 
Points distributed along the main parabolic arc 
are due to interference between a highly scat- 
tered component at 6 and the components aris- 
ing from near the pulsar brightness centroid at 
= 0. In this case there is a differential delay 
r = 2 D c g/2c and differential Doppler frequencjo 
/d = V e ff • 8/ A. Here D D g and V c g are the effec- 
tive distance and velocity vector of the scattering 
screen as defined in 21 and A is the observed wave- 
length. If the scattered image is essentially linear, 
through the origin at angle a with respect to V e ff, 
we can eliminate 9 and obtain 

r = a/n , where a = D cS A 2 /(2cVjf cos 2 a) (1) 

as the equation of the main parabola. Note that 
a oc A 2 regardless of the wavelength scaling of 
scattering angle 9, as confirmed experimentally by 
Hill et al.1 (|2003ll . 

To understand the interferometric observations 
consider the pulsar as a point source at a dis- 



1 Henceforth these terms will be abbreviated to delay and 
Doppler. 
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Differential Doppler frequency f (mHz) Differential Doppler frequency f Q (mHz) 

Fig. 1. — The amplitude (left) and phase (right) of the seconadary cross spectrum C, defined in Table Q] for 
the 8 MHz wide 314.5 MHz sub-band on the AO-GBT baseline. Amplitude scale is log 10 and phase is in 
degrees For display purposes only, both diagrams have been smoothed over 5 pixels in delay to reduce the 
noise, giving plotted resolution of 0.44 mHz in Doppler frequency and 3.1 fis in delay. The black line is a 
parabola with curvature 0.56 s 3 on both panels. There are bright arclets in amplitude that extend along the 
ridge of this main parabola to delays > 1 ms. The phase appears much smoother, as discussed in jj3.ll 



tance D p from Earth, whose radiation is scat- 
tered by an inhomogeneous region of ionized in- 
terstellar medium at a distance D s from Earth. 
In this thin-screen approximation the phasor due 
to radiation scattered at a point x.,- = D s 9j in 
the screen suffers a phase delay 4>j due to the 
scattering plasma and is received with a total 
phase delay <F, (r) — (\> } + fc(xj — j3v) 2 /2(3D S where 
j3 = 1 — D s /Dp and r is the transverse position 
of the receiving antenna. The received electric 
field is the summation of all such scattered compo- 
nents, either expressed by the Fresnel diffraction 
integral (see lCordes et al.l 120061 ) or approximated 
by a summation over a ll stationary phase points 
see IWalker et al1l2004h . in strong scattering. 

Here we use the stationary phase approxima- 
tion to obtain the correlation of the fields between 
antennas at locations — b/2 and b/2. This defines 
the visibility whose Fourier transform in time and 



frequency can be written as: 

V(r,fn,b) = 5> Xp [;(%(-b/2) - %(b/2))] 

X MjMfe<J(/D - fu,jk)S(T - Tj k ). (2) 



where: 



Tjk 



{dj - e k ) • v eff , 



<Pj _ 9k 

2ixv 2ttv 



(3a) 
(3b) 



The S functions in equation ^ determine how 
each particular point on the secondary spectrum 
is related to a pair of points in the screen. With 
finite bandwidth and observing time the S func- 
tions should be replaced by finite narrow sine func- 
tions. The equation sums over all possible pairs of 
stationary phase points in the screen, with Hj,fJ>k 
as the magnifications de termined by phase cur- 
vature of each point (see Walker et al. 2004). In 
the full Fresnel formulation, the summation be- 
comes a double two-dimensional integration over 
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Term 


Symbol 


Remarks 


Correlator products (2nd order in electric field) 


Dynamic spectrum 


Kf,t) = (\E(f,t)n 


the power of the electric field as a function 
of frequency and time 


Dynamic cross spectrum 


V(f,t,b) = (Ei(f,t)EZ(f,t)) 


the visibilities measured on a baseline b as a 




function of frequency and time 


Intermediate quantities 


(2nd order in electric field) 






/(t,/d) = FT [J(/,t)] 


the Fourier transform of the dynamic spectrum 




V{r,f D ,h) = FT[V(/,t,b)] 


the Fourier transform of the dynamic cross spectrum 


Derived quantities (4th 


order in electric field) 




Secondary spectrum 


A(r,fn) = \I(r,fn)\ 2 


the squarcd-modulus of the Fourier transform of the dynamic 
spectrum 


Secondary cross spectrum 


C(T,/ D ,b) = 


the (complex) product of V(t, /.D,b) and its 




V{r,f D ,h)V{-T,-f D ,h) 


corresponding quantity measured at the point ( — r, —fo) 



Table 1: Definitions of various products derived from the measured intensities and visibilities used throughout 
this text. In the definitions for the first two quantities, E is the electric field as measured by the antenna 
receiver and the angle brackets denote a time average over the measurement interval much longer than the 
inverse bandwidth. 



the screen with all ji = 1. The plasma delay term 



in square brackets in equation (|3b[) involving <fij is 
unimportant compared to the first term when con- 
structing secondary spectra in strong scintillation, 
but in any case this term cancels in the astrometry 
discussed below. 

In either formulation equations (|3a|) and (|3b[) 
connect the points on the scattering screen to r, /d 
in the secondary spectrum and are independent of 
r. Note that /d can be understood as the differen- 
tial Doppler shift between tw o waves scattered at 
differing points on the screen ( Cordes et al. 20061 ) 
and that the delay r is the differential group delay 
evaluated at the center of the observing band. 

Each arclet arises due to the interference of ra- 
diation from a given fixed point on the scatter- 
ing screen with scattered radiation from 8^0. 
Examination of the delay in equation (|3b[) shows 
that for a strong contribution from fixed 6j the 
greatest delay (i.e., arclet apex) will be where Ok 
goes through zero. The particular shape of each 
of the arclets reflects the scattered brightness dis- 
tribution near core of the image. Furthermore the 
fact that the SS amplitude is bright over a narrow 
range in delay shows that the points comprising 
the bright emission centroid region are extended 
along a line rather than in a circula r halo (see 



Cordes et al.ll2006l; IWalker et al.ll2004l) 



The velocity of the scintillation pattern across 
the line of sight is V ISS = (1 - /3)V p + /3V e - V s , 



where V s , V 

screen, Earth and pul sar respectively (see e.g 



and V p are the velocities of the 



Cordes fc Rickettlll998h . The scintillation veloc- 
ity is dominated by the pulsar, whose proper 



motion ( Lvne et all Il982t ) combined with a dis- 
tance of 643 pc (estimated from the dispersion 
measure using the NE2001 Galactic electr on dis- 
tribution model of ICordes fc La"ziol (l2002t n gives 



p. a 



6 1+ 15 

o.i_ 15 



km 



p,6 



km 



note that these velocity errors do not include the 
substantial uncertainties in the pulsar distance. 

3.1. Astrometric Imaging 

In single dish observations information on the 
scattered image is limited to the power in the sec- 
ondary spectrum, obtained from I(t, /d). How- 
ever, the addition of interferometric observations 
permits high-precision astrometry to be performed 
on each component of the secondary cross spec- 
trum, which isolates the wavefield due to a single 
pair of interfering waves. By applying astrometry 
to each such pair individually, we can recover an 
apparent image of the scattered pulsar radiation. 

For a particular pixel from V(t, /d, b) equation 
([2]) sums over all pairs of points that satisfy equa- 
tions ((3]). However, we now consider the case that 
only one pair of points satisfies the latter condi- 
tions, a condition that is valid for a linear or suf- 
ficiently anisotropic source. Then the phase of V 
is$ ilfc = $j(-b/2)-$ fc (b/2)is 



kD s 

w 



/3b 



(fli+e k )] (4) 



which has terms that are symmetric and anti- 
symmetric in scattered points j, k. We cancel the 
anti-symmetric terms by summing the phase of 
V at each point (/d,t) and its conjugate point 
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at (-/d,-i"). This eliminates the random screen 
phase leaving the symmetric part of the phase as 

Vj,fe(b) - <!>,., + <l>,,, - • ■(), + 0, i , (5) 

which neatly encodes the average position of the 
two scattering points projected parallel to the 
baseline b. As we show below, the scattering is 
indeed highly anisotropic and the assumption of 
a single pair of interfering waves for each pixel is 
well justified. 

In practice the anti-symmetric phases are can- 
celled by first forming the secondary cross spec- 
trum, defined as 

C(/ D , t, b) - V(f D ,r, b)V(-f D , -r, b) , (6) 

for all non-negative values of r. Subsequently, the 
resultant complex phasors are averaged over small 
regions (3 pixels in Doppler by 5 in delay) in order 
to reduce the effect of noise before the argument 
is determined: 



V(b) = arg«C(/ D ,T,b)» 



(7) 



Here () denotes the 3 by 5 averaging which prop- 
erly weights the complex products before the 
phase is computed. The result is the geometric 
phase of the secondary cross spectrum sampled 
with a resolution of 0.63 fis in delay and 0.44 mHz 
in Doppler. As with any phase measurement the 
result is modulo 2ir. 

The right panel of Figure [1] shows tp, the phase 
of the secondary cross spectrum, for the baseline 
from GBT to AO. In contrast to the sharply de- 
fined arclets in the amplitude plot, the phase is 
relatively smooth, appearing continuous between 
the arclets, as might be expected if the fine struc- 
ture in amplitude comes from a continuous image 
with fine structure in brightness. Consequently we 
can use astrometry from tp to actually map out the 
scattered brightness distribution with remarkable 
precision. Furthermore we can measure the posi- 
tion angle of the axis of elongation and estimate 
both the distance to the region of scattering and 
its velocity. 

We can use equation ([5]) to find (Gj + Ok) for 
any point in the secondary cross spectrum if it 
is dominated by a single pair (0j,8k) and there 
are at least two baselines with sufficient signal- 
to- noise ratio (S/N). However if we can identify 



the apex of an arclet we know 6j or 6k is zero, 
so we have a unique solution for the remaining 
angle. The amplitude plot in Figure [1] shows that 
the apexes of the identifiable arclets lie near the 
parabolic curve shown as a line. Hence we sampled 
the secondary cross spectrum along that curve to 
approximate a set of apex positions even where 
the arclets are too densely spaced to be identified 
individually. We defined the main parabolic arc by 
the curvature a = 5.577 x 10 16 // 2 (s 3 ) for center 
frequency / Hz in each of our four sub-bands, as 
shown by the black line in Figure [TJ 

We selected data with a high enough S/N to 
estimate 6 as follows. For each of the six baselines 
we found the phasor (C(/d, t, b)) and defined S as 
the S/N relative to the root mean square value of 
| (C) | in a low power (noise dominated) region near 
the point being examined. Samples with delays 
within ±0.6 /is of the main parabola were selected 
if S > 4 on the GBT- AO baseline and if S > 3 on 
the WB-AO baseline. 

For each sample of (C) we first formed a "dirty 
image" from the 6 baselines in the same sense as 



for tr aditional synthesis imaging (jTavlor. Carilli fc Perley 
1999). At the dirty image maximum a two- 



dimensional Gaussian was fit to determine a po- 
sition, and both S and the shape of the Gaussian 
were used to assign positional uncertainties. How- 
ever, since C is a product of two visibilities this 
method gave too much relative weight to the high 
S/N baselines and led to apparent positional er- 
rors much smaller than the scatter of the points 
from neighbouring samples along the main arc. 
In addition there are lobe ambiguities in the syn- 
thesized beam, since two of four antennas (JB 
and WB) are relatively close to each other and 
so probe nearly the same spatial frequencies. The 
result is a dirty beam that is nearly doubly peri- 
odic with lattice vectors a! = (59.0,-61.9) mas 
and a 2 = (34.2, 17.8) mas at 322.5 MHz (scaling 
linearly with wavelength). 

After some investigation we adopted a simpler 
astrometric method based on a weighted linear 
least squares fit of the observed phases, tp, using 
equation © (with one of the angles set to zero). 
Under conditions of high S/N one can show that 
the phase error is cx 1/ ^/S and so we used weights 
oc VS". Because the phase is observed modulo 2tt 
there are ambiguities in position, corresponding to 
the lobe ambiguities. Consequently, we used a first 
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Relative Right Ascension (mas) 



Fig. 2. — Map of the astrometric positions of the 
scattered radiation from samples along the main 
parabola in the sub-band at 322.5 MHz plotted as 
error bars in colors of green, orange, red & dark 
red, for amplitude, |C|, increasing by factors of 10. 
Right ascension and declination values are relative 
to the centroid position as determined by VLBI 
self-calibration. Note the increase of astrometric 
errors with decreasing amplitude \C\. The blue 
straight line is an unweighted fit to these points. 
The black error bars are from the apexes of iden- 
tified arclets and the black line is a fit to these 
points. Samples from the 1 ms feature (negative 
Doppler and delay > 0.9 ms) are shown by cyan 
points for two possible lobe choices. As discussed 
in 2] the correct lobe position is at the lower right. 
This is overplotted with black error bars from the 
apexes of individual identified arclets near 1 ms. 



guess position as an initial model and wrapped 
the phase (e.g., on the transatlantic baselines) to 
be within ±7r of that solution before applying the 
least squares fit. For the points from the main arc 
our first guess was at the origin and the results are 
plotted as error bars (in a pale color) in Figure [2] 
Points from the main arc form a continuous elon- 
gated image close enough to the origin that the 
phase ambiguity is not a problem for them. 

An empirically determined line passing through 
arclet apexes was used to sample the points in 
the SS associated with the 1 ms feature. The 
astrometry of these points, however, forms a dis- 
tinct group offset from the origin in celestial co- 
ordinates, making it difficult to resolve the phase 
ambiguities. Accordingly, we used first guess posi- 
tions from the two closest synthesized beam lobes 
which are roughly equidistant from the origin. The 
cyan points on the left show the fitted positions 
starting from the central lobe; the cyan points 
on the right show them starting from a position 
lobe-shifted by the lattice vector —a.2, which cor- 
responds to the addition of 27r to the phase on the 
transatlantic baselines. In S|4]we discuss the reso- 
lution of this ambiguity in favor of the right hand 
position. 

The major conclusion from Figure [2] is that the 
scattered image from the main arc is remarkably 
extended along a straight line with very few points 
more than 2cr from that line. As noted above 
highly anisotropic scattering was already inferred 
by the emptiness of the secondary spectrum (par- 
ticularly the absense of power at fr> ~ 0, r > 0) 
and the many discrete arclets. The apparent flar- 
ing of the points away from the origin is due to 
the increasing size of the error bars because the 
SS amplitude decreases with angle of scattering. 
We define the straight line fitted to the points as 
the "scattering axis" . Note that the plot is not a 
proper image of the scattered brightness, since the 
amplitude |C| must be transformed by a Jacobian 
into the brightness. Hence the axial ratio is not 
easily estimated and we postpone the discussion 
of it to a later section. 

3.2. Arclet Apex Astrometry 

The astrometric positions discussed in the pre- 
ceding section and plotted in Figure[2]are based on 
the assumption that samples along the main arc 
correspond to the apexes of an underlying arclet. 
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However this is not always true, as is obvious in the 
1 ms feature. So we have made an effort to identify 
all the discrete arclets, estimate the apex of each 
in an optimal way, and then map the phases at 
each apex to an angular position. At 322.5 MHz 
we found 62 well defined arclets (with apex delays 
t > 0.1 ms) of which 7 are in the 1 ms feature. 
At smaller delays the arclets are too crowded to 
identify separately. 

We formed a template by averaging all the ar- 
clet amplitudes, centered by their apex, and we 
fitted a straight line in phase versus /d to the 
template to obtain the phase gradient versus /d in 
the template. We then fitted the amplitude and 
phase templates to each arclet using a weighted 
least squares fit. Hence we obtained estimates for 
the delay, Doppler, phase, and amplitude at the 
apex and errors in all four parameters. We then 
estimated 6 and its error at each apex and re- 
peated the procedure for all 4 sub-bands. 

The astrometric positions and errors of the 62 
arclets are plotted in black in Figure [2] The error 
bars are smaller than from the main arc because 
the phases were determined from a weighted av- 
erage of the phases over the entire arclet. In ad- 
dition the scatter should be smaller because they 
are known to correspond to interference with the 
bright centroid component at the origin. How- 
ever, the lobe uncertainty in the location of the 
7 arclets from the 1 ms feature remains. As we 
discuss below, we resolve the ambiguity in favor of 
the lower right lobe position marked by the black 
error bars derived from the 7 apexes. The apexes 
near the main arc have positions extending out to 
20 mas from the origin and show good linear align- 
ment with the inner points from the main arc. A 
straight line fitted to the inner main arc has a po- 
sition angle of — 25.2±0.5° (east of the declination 
axisjl. We use this axis to define 9\\ and define 9± 
by a 90° clockwise rotation from it. 

We have used three techniques to resolve the 
lobe ambiguity. The first relies on the frequency 
dependence of the positions obtained from each 
arclet apex. As discussed below in ^3 . 31 we find 
that the angular positions of the scattered fea- 
tures that cause the arclets are essentially inde- 
pendent of frequency. Thus we can use the fact 



that the correct lobe choice should give a position 
that is independent of frequency. We computed 
the 2 dimensional mean for the 1 ms apexes from 
a weighted average position of its imaged points 
separately for each of the four sub-bands. For the 
lobe position at a negative RA offset in Figure [2] 
the mean position was independent of frequency 
within the errors of estimation (~ 0.1 mas). All 
other lobe choices show systematic frequency de- 
pendence in their positions in agreement with that 
expected for a lobe error. 

Confirmation comes from two further tests. 
The second technique relies on equation (|3bj) relat- 
ing delay to the angular offset of an arclet apex and 
is discussed in H4.ll Similarly, the third technique 
relies on equation (|3a[) for the Doppler frequency 
and is discussed in £14.21 



3.3. Frequency Scaling 

Hill et alJ (|2005h found the wavelength scaling 



of the apex Doppler frequencies of the four arclets 
which they identified in pulsar B0834+06 to scale 
oc A -1 . Hence using equation (|3a|) they concluded 
that the scattering angle responsible for the arclet 
was independent of wavelength. We have done a 
similar analysis on the apex Doppler frequencies of 
all the identified arclets in our four sub-bands. In 
addition we analyzed the wavelength scaling of the 
apex delay in the four sub-bands, which is oc 9 2 as 
in equation (|3b|) . We combined the fits to estimate 
a single scaling parameter 7 where 9 oc A 7 . We 
estimated 7 separately for three groups of apexes 
and tabulate the values in Table O 

Table 2: Wavelength-scaling exponents, 7, in the 
scattering angles (9 oc A 7 ) estimated from apex 
positions for three groups of arclets. 
Arclet apex group 7 



r ~ 1 ms /d < 
0.1 < t < 0.4 ms 
t > 0.4 ms f B > 



0.062 ±0.006 
0.01 ±0.01 
0.019 ±0.004 



2 Here and elsewhere in this paper, reported uncertainties 
are 68% confidence intervals. 



The results show that the scattering responsi- 
ble for each arclet originates at a location that is 
essentially independent of wavelength across the 
10% range spanned by the four sub-bands. This is 
entirely incompatible with the value 7 = 2, to be 
expected if the individual arclets came from sepa- 
rate ray paths caused by refraction due to differing 
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gradients in the column density of electrons. The 
conclusion applies to all arclets in the SS, with 
the possible exception for the 1 ms feature whose 
position scales weakly A 06 . Such a weak scaling 
might apply if the arclet were caused by a lens-like 
concentration of pla sma, such as that invoked by 
(jRomani et al.l[l987h to explain ESEs, but with a 
transverse extent very much smaller than than its 
distance from the direct path. 

4. Scattering Model 

We now develop a model for the distance and 
velocity for the various features in the scattered 
image. 




Fig. 3. — 9 2 versus delay. Points from the main 
arc in Figure [2] are shown as blue x (/b > 0) and 
o (/d < 0) ; the apexes of the identified arclets are 
plotted in black, except for those from the 1 ms 
feature which are in red. The straight line models 
(9 2 = mr; see equation©) were fitted separately 
to various subsets of the points to estimate their 
distances. 

4.1. Scattering Distance Estimation 

Equation (|3b|) relates delay r to the angular 
offset 9 for an arclet apex, for which the second 
angular component is at the origin. Thus: 

r = (D cS /2c)9 2 (8) 



where 

D cS = D p (l - 13)/ P (9) 

In Figure [3] we plot on logarithmic scales 9 2 from 
the astrometry against delay with the various 
groups of data from Figure [2j Note that there is a 
bias to 9 2 due to the uncertainties in right ascen- 
sion and declination. Since we have estimates of 
those uncertainties we subtracted the bias before 
fitting. 

Two separate straight lines through the origin 
were fitted to the (blue) points sampled along the 
main parabola and from the 1 ms feature. The 
apex points were also fitted separately for three 
regions: for apex delays less than 0.4 ms, delays 
bigger than 0.4 ms with positive Doppler and with 
negative Doppler (i.e., the 1 ms feature). The fits 
were weighted proportional to the reciprocal of the 
standard deviation in 9 2 . The derived effective 
distances are given in Table [3] with their formal 
la errors for all 4 sub-bands. 

The simplest model puts all of the scattering 
in a single relatively thin region at the same ef- 
fective distance. Comparing D B g for the differing 
regions of the SS, we find D s in the range 1100 to 
1400 pc. Since we found no significant scaling of 
the apex 9 with frequency (see $3j3j), we combine 
data from all four sub-bands to create a combined 
estimate shown as a fifth row in the table. We 
note that there is a possible contribution to in- 
consistencies among the sub-bands from a partial 
narrowing of the two outer sub-bands at AO by a 
front-end bandpass filter. 

The average of the estimates in column 2 from 
delays less than 0.4 ms from both the main arc and 
from the arclet apexes is 1171 ± 23 pc. In column 
4 the 1 ms cluster gives 1301 ± 40 pc, which over- 
laps at ±2<T with the value 1121 ± 59 pc obtained 
from the 1 ms apexes. This small discrepancy ap- 
pears in Figure [2] as an offset in the centroid of 
the blue and black points. Since the arclet apexes 
have been specifically identified so that one of the 
6 parameters can confidently be set to zero and 
since D c g for the 1 ms apexes is within la of those 
from the main arc, we conclude that the scatter- 
ing region responsible the 1 ms feature is located 
at the same distance as that causing the main arc. 
Without the lobe-shift a fit for D c g gives values 
from 1700 to 2600 pc. The large discrepancy in 
this distance provides the second confirmation of 
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Table 3 

Effective Distances and Velocities from the VLBI astrometry 



Main Parabolic Arc 


Frequency 




D e{t (pc) 


y. ff (km s- 1 ) 


MHz 


Inner Arc 


r ^ 1 ms 


K tf || a 


V ef fX b 


314.5 
322.5 
330.5 
338.5 


1197 ± 23 
1175 ± 24 
1195 ± 29 
1133 ± 18 


1350 ± 65 
1272 ± 29 
1294 ± 20 
1287 ± 24 


303.3 ± 3.1 
300.7 ± 3.3 
303.3 ± 3.7 
300.0 ± 2.5 


-131.9 ± 6.4 
-151.5 ± 11 
-161.5 ± 16 
-133.8 ± 5.0 


326. 5 C 


1175 ± 29 


1301 ± 40 


301.8 ± 3.2 


-144.7 ± 12 


Apexes of Identified Arclets 


Frequency 




D eti (pc) 


V b{{ (km s _1 ) 


MHz 


0.1 < r < 0.4 ms 


r > 0.4 ms /d > T ~ 1 ms 


Veff|| a 




314.5 
322.5 
330.5 
338.5 


1165 ± 19 
1205 ± 19 
1162 ± 17 
1139 ± 15 


1370 ± 43 1064 ± 63 
1419 ± 67 1132 ± 24 
1406 ± 39 1180 ± 110 
1315 ± 67 1107 ± 17 


305.9 ± 2.9 
313.4 ± 2.5 
306.8 ± 2.3 
305.4 ± 2.4 


-136.1 ± 4.4 
-152.9 ± 3.7 
-158.4 ± 11 
-134.9 ± 1.8 


326.5" 


1168 ± 23 


1378 ± 60 1121 ± 59 


307.9 ± 3.3 


-145.6 ± 8.6 


a Estimates 


of Veffy were derived from the main arc, since they lie along 


the scattering axis 





b Estimates of V e {f± were derived from the 1 ms points, by assuming that they have the same V e ff|| as for 
the main arc 



c Averages from all four sub-bands with errors that include the uncertainties in each sub-band. 



the lobe-shift for the 1 ms points. 

A possibly significant difference is seen between 
the apexes with positive Doppler and delays > 
0.4 ms in column 3, which give a significantly 
larger D g = 1378 ± 60 pc. Whereas all of the 
sub-bands give a larger distance there are only a 
few identified arclets in each sub-band and they 
are of relatively low S/N in all 4 sub-bands so 
that accurate correction for the noise bias to 2 
becomes more critical. Thus the evidence is rela- 
tively weak for a different distance for these points 
at large delay and positive Doppler. 

We now take D e ff = 1171 pc, combined with 
the pulsar distance of 640 pc in equation © 
to obtain the fractional scattering distance (3 — 
0.353 ± 0.005. This is consistent within 2a of the 
value of 0.29 ± 0.04 obtained bv lHill et al.1 (|2005l ) 
using the same pulsar distance but without the ad- 
vantage of knowing the full geometry of the scat- 
tering screen. However, its true uncertainty is 
dominated by uncertainty in the pulsar distance 
which is estimated from the dispersion measure. 
Dispersion measure distances are notoriously un- 



reliable, especially for ne arby pulsars and thu s are 
only accurate to ~ 40% (jBrisken et al.ll2002f ) ; this 
gives a screen distance range of 250 to 580 pc from 
the Earth. 

4.2. Velocity Estimation 

An independent analysis is provided by equa- 
tion (|3a|) . relating the Doppler frequency for the 
arclet apexes to their astrometric position. At 
each apex the Doppler frequency can be written 
as: 

/d = 0-V cS /X 

= (0\\V cSll +e ± v cS1 _)/x. (io) 



where 



v„ 



Viss//3« V p (l-j0)//J. 



(11) 



The final approximation ignores the velocities of 
the Earth and of the screen relative to that of the 
pulsar. Then the ratio V cr f/Z3 c ff estimates the pul- 
sar proper motion independent of its distance and 
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Doppler (mHz) 

Fig. 4. — 9ii plotted against the Doppler /rj for 
the sub-band at 322.5 MHz. Points plotted with 
cyan x marks are sampled from the inner main 
arc in Figure [2] and from a line through the 1 ms 
feature. Black points are from apexes of identified 
arclets, except 1 ms apexes which are red. 

of j3. However, here we do not make that approx- 
imation but use the measured proper motion to 
estimate the velocity of the scattering screen. 

In Figure [3] we plot 6n against /b which would 
give a straight line whose slope gives an estimate 
of A/14ff || , providing that the value of 9±V c s± were 
the same or negligible for all points. As shown in 
Figure [2] the points fitted by the straight line have 
small 9± scattered about zero. Hence we were able 
to estimate T4ff|| from the various groups of points 
and tabulate the results in Table [3] excluding all 
the points with delays > 0.4 ms. 

If we combine inner apex fits with values from 
the main arc we obtain V c g\\ = 305 ± 3 km s _1 . 
We now assume the same V e g\\ to apply to the 
1 ms apexes and since these points have a signifi- 
cant measured 9j_ we can fit for the perpendicular 
velocity also, obtaining V c g± = —145 ± 9 km s _1 . 

In summary the scattering model provides three 
observable parameters which we list in Table [4] 
We can express the velocities in RA, Dec and re- 
late them to the measured pulsar proper motion 



velocity V Pta = 6.1 km s _1 , V Pl s — 156 km s _1 . 
The derived model parameters are also listed in 
the table. Equation (jTTJ) relates the measured ef- 
fective velocity to the scintillation velocity Viss 
which also depends on the transverse projected 
velocity of the Earth and the scattering screen. 
Using the pulsar velocity cited above gives an es- 
timate of the ISM screen velocity relative to the 
Sun (~ 16 km s _1 ), which is comparable in mag- 
nitude to other measured interstellar velocities. 

The ISM velocity analysis provides the third 
confirmation of the lobe-shift applied to the 1 ms 
astrometry, since in the alternate position V c f[± 
is reversed and could only be reconciled with the 
pulsar velocity if there were an implausibly high 
velocity (> 200 km s^ 1 ) for the scattering screen 
relative to the Sun. 

Table 4: Model Parameters for Distance and Ve- 
locity, assuming (3 — 0.353. Note that the first 
five quantities are measured; the others are calcu- 
lated from these assuming the pulsar distance and 



velocity as cited in 


the text. 








1171 ±23 


pc 




Ken 


305 ±3 


km 


s- 1 


v cS± 


-145 ±9 


km 


s- 1 


| Scattering axis 


-25.2 ±0.5 


deg 


. east of north 


_L Scattering axis 


-115. 2± 0.5 


deg 


. east of north 




415 ±5 


pc 






-16 ± 10 


km 


s- 1 


v s± 


0.5 ± 10 


km 


s- 1 


a 


27 ±2 


deg 





"Assumes D p = 640 pc; the error in D B due to the uncertainty 

in Dp is much larger ~ 40% 
^Including errors from uncertainty in pulsar distance and 

proper motion 

4.3. Estimating the Image Center 

In the discussion so far we have assumed that 
the astrometry is correctly centered on the emis- 
sion centroid. However, since the scattered im- 
age is highly anisotropic, the centering done by 
self-calibration in the primary analysis would have 
been more precise in the perpendicular direction 
than in the parallel direction. Close examination 
of the initial versions of the plots as in Figures [3] 
and [4] revealed small but significant asymmetries, 
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Fig. 5. — Cyan error bars are the VLBI astrom- 
etry of the apexes of the identified arclets super- 
imposed from all 4 sub-bands. Red error boxes 
show the much finer precision of the back-mapped 
astrometry, as explained in i j4.4l and clearly de- 
lineate the highly linear scattering from both the 
main arc and the 1 ms feature (at lower right). 

which allowed us to improve the centering in the 
parallel direction. 

In the Doppler plot of Figure [4] a straight line 
fit did not pass exactly through the origin and this 
fitted offset in 6u was found to be on the order of 
0.5 mas differing slightly between sub-bands. In 
the delay plot of Figure [3] the 9 2 values at a given 
delay with positive Doppler frequencies were sys- 
tematically shifted relative to those with negative 
Doppler frequencies. We were able to minimize the 
X 2 for the straight line by fitting for a positional 
offset in the parallel direction. We found that in 
each sub-band the optimum parallel offsets were 
consistent between the Doppler and delay estima- 
tion methods within their errors. An average of 
these two offsets was consequently applied to the 
astrometry in each sub-band before the analysis 
described in the two foregoing sections. 



4.4. Astrometry Mapped Back from the 
Secondary Spectrum 

The synthesized beam size of our VLBI array 
(its basic angular resolution) was about 35 mas at 
327 MHz. In the astrometry shown in Figure[2]we 
obtain positional accuracies as small as ±0.5 mas 
near the center of the image due to the high signal 
to noise ratio and much larger errors away from 
the center. However we can do even better from 
the apex positions of the arclets in the SS using 
equations ijHJ) and (Ti"0|) . which we will refer to as 
"back-mapping". This is possible since we used 
the VLBI astrometry to obtain accurate values for 
the model parameters D e g , V e g, and a. The equa- 
tions give the magnitude of 8 from the delay r and 
the component 8y parallel to Viss from /d- There 
is an inherent ambiguity in such a mapping since 
it does not determine the sign of the angular co- 
ordinate orthogonal to the velocity. However, we 
can use the VLBI image to resolve that ambigu- 
ity and show the positions of each apex mapped 
back into in RA, Dec in Figure [5] overplotted for 
the 4 sub-bands. The mapped astrometry points 
are plotted as red error boxes on a background of 
the error bars of the VBLI positions. The apex 
points mapped back in this fashion extend along 
the same axis as the VLBI positions but with a 
narrower spread in 8± due to their even smaller 
errors. 

With errors St and Sfx> in the apex positions we 
get positional errors S9 ~ cSt/(D c s8) and 88y ~ 
Sfo^/Veff. With St ~ l we obtain 89 ~ 8^ asl 
i.e., 0.05 mas at # mas = 20 mas, which improves on 
the 0.5 mas from astrometry, particularly at large 
scattering angles where the S/N worsens; and with 
Sf-D ~ 1 mHz the resolution S8y ~ 0.6 mas. The 
resulting confidence region is shown by the red box 
which is narrow in radial 8 coordinate and wide in 
the transverse direction. 

The map emphasizes the remarkable anisotropy 
of the scattering disk and also shows evidence that 
the 1 ms feature is also highly extended along a 
roughly parallel direction. This is an important 
aspect that will constrain any physical model for 
the scattering. 

4.5. Scattered Brightness Function 

The scattered brightness distribution can be re- 
covered from the secondary spectrum if the intcn- 
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Fig. 6. — Le/f;Scattered brightness against 9u obtained for points along the main arc via the back-mapped 
astrometry in £ 14.41 averaged from all four sub-bands (blue). The individual peaks are as narrow as 0.1 mas 
as shown in the expanded view in the right panel. The three overplotted theoretical curves (red) are for a 
one-dimensional Kolmogorov model. The middle curve was fitted to the observations over the range shown 
in green. The other two curves have the same total flux density for the pulsar but are wider and narrower 
as discussed in the text. 



sity scintillat ions are weak, or if they are highly 
anisotropic ( Cordes et al. 2006). In this case the 
scintillations are highly anisotropic and we can 
make one-dimensional, strip integrated brightness 
distributions, Bi(9u), by strip integrating the two 
dimensional distribution over 9\ . Then we have 



Bi(9 1 )Bifo)=A(r ) f D )\J\ 



(12) 



where 



0j)D cS /2c, and / D = (0 1 - 



92)V c ff cos a/ A and the Jacobian |J| = ,/b-Deff/c. 
Then by sampling A(t, /d) along the main arc 
where 62 = 0, we can estimate B\{9). This es- 
timate is plotted in Figure 6 without the 1 ms 
feature, which cannot be represented on this plot 
since it is not on the main parabola. 

If the scattering comes from homogeneous an- 
isotropic Kolmogorov turbulence then the strip 
integrated brightness distribution is the Fourier 
transform of the one-dimensional correlation 
function of the electric field, Bi(9\\) = FT 

[exp(— (sy/so) 5 ^ 3 )]- We have plotted three such 
models over the observations in Figure [H each 
with a different value of sq, the coherence scale of 
the electric field. The middle model fits the data 
near the origin (so ~ 10 4 km) and the other two 
give a rough estimate of the range of sq necessary 
to match the measurements. They correspond to 



changing the root mean square (RMS) electron 
density by a factor of 2 (wider curve) or 0.5 (nar- 
rower curve). One can see in the expanded view 
on the right that near 9\\ = -10 mas there is a 
change in the RMS electron density by a factor of 
4 in a very small distance. 

The finest structure has an angular scale of 
about 0.1 mas, which corresponds to 0.05 AU. This 
plot can also be made using the VLBI astromet- 
ric 6*|| . In this case the large scale structure is the 
same but the small scale structures seen clearly in 
Figure [S] are smoothed out by the larger error bars 
in the VLBI astrometry. 

4.6. Axial Ratio 

While we have clear evidence for highly aniso- 
tropic scattering, it is difficult to estimate the axial 
ratio R, which is usually defined by a contour at, 
say, half power in the two-dimensional scattered 
brightness distribution. Estimates of the width in 
#11 can be made from the one-dimensional bright- 
ness in Figure [5] However, the width at the same 
level in 9± is harder to estimate, since we have no 
information on brightness versus 9±. As an alter- 
native, we can define an apparent axial ratio i? ap 
from the scatter of the points in Figure [5] indepen- 
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dent of their brightness. Here we define R ap as the 
ratio of the RMS width in 6±_ to that in 0y. The 
observed scatter in 9± will of course be broadened 
by the astrometric errors, so that only an upper 
bound can be found on 9± due to interstellar scat- 
tering. 

We apply these ideas to Figure [71 which shows 
the astrometry from the main arc as a scatter plot 
of 9± against 9u superimposed from all four sub- 
bands. The black error bars show 9± averaged into 
0.5 mas bins in 9\\ where the length of the bar is 
the standard deviation in each average. Taken as a 
group the black points from the inner ±5 mas are 
consistent with a Gaussian distribution with zero 
mean and 0.38 mas RMS; there are no points more 
than ±2cr from zero and their error bars, which 
have not been corrected for the measurement er- 
ror, have a mean value of 0.3 mas. Thus the true 
scatter in 9± is smaller than these errors, implying 
an intrinsic perpendicular RMS width less than 
~ 0.3 mas. Taking the half width at half power in 
9u to be 3 mas from Figure [6] gives a lower bound 
on the axial ratio R ;> 10. Similar estimates can 
be made from the ratios of the parallel to the per- 
pendicular RMS widths of the apex astrometry 
from both the VLBI and back-mapping methods. 
These yield R ap J> 27 and i? ap J> 20, respectively. 

The scattering axis is neither parallel nor per- 
pendicular to the Galactic plane. However, ac- 
cording to the Viriginia Tech Spectral-Line Sur- 
vejH the pulsar lies ~ 1° from a 5° long bright 
Ha filament at a position angle within about 10° 
of our scattering axis. Although there is a rough 
agreement in position angle we have no other ev- 
idence to support an association with our 30 mas 
long filament. 

5. Physical Models of the Scattering 

While our main focus is the presentation of the 
remarkable results found from the observations, 
we now briefly consider what physical structures 
could be responsible for the scattering. Following 
previous analyses of the ISS phenomena, we as- 
sume that the basic cause is scattering by random 
(presumably turbulent) structures in the plasma 
density. The scale of such turbulence must extend 
down to well below the diffractive scale, ~ \/2it9\\ . 

3 See http://www.phys.vt.edu/halpha/ for more informa- 
tion about this survey. 



To obtain scattering angles of 25 mas requires mi- 
crostructure of the order of 1000 km. However, 
the turbulence is inhomogeneous over scales of 0.1 
to 10 AU as mentioned above and we must also ac- 
count for the highly elongated shape of the main 
scattered image. 
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Fig. 7. — Astrometry referred to axes parallel and 
perpendicular to the main scattering axis. Pale 
color superimposes scatter plots from the main arc 
in all four sub-bands. The black points show the 
average 9± from bins 0.5 mas wide in 9\\ with ver- 
tical bars giving the error in the average. The bars 
have not been corrected for the astrometric errors. 
Note that only the inner portion of the primary 
scattering disk where high S/N points are located 
is included in this plot. 

First we ask whether the peaks in our scattered 
image, due to the arclets, could be the equiva- 
lent of "speckles". Discrete arclets are visible in 
the simulated secondary spectrum in Figure 11 of 
Cordes et al. (2006), who computed the SS for 
a realization of a layer with a Kolmogorov spec- 
trum and an axial ratio 4:1 parallel to the simu- 
lated velocity. That figure shows very many fine 
arclets with apexes centered on a pronounced for- 
ward parabola. In independent realizations the ar- 
clet apexes appear at random locations, showing 
that they are the equivalent of speckles in a snap- 
shot scattered image through a turbulent scatter- 
ing layer. Thus they do not map to deterministic 
structures in the layer. 

In contrast the discrete nature of the arclets 
and their sparser distribution in our observation 
suggests deflection by discrete structures, possibly 
similar to those responsible for extreme scattering 
events (ESEs). The two basic models proposed 
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for ESEs have been an e nhancement in n R t hat 
causes a diverging "lens" (jRomani et al.lll987j ) or 
an enhancement in plasma turbulenc e that causes 



extra scattering (jFiedler et al.l 119871 ). The result 



that the angle of arrival from the arclets is essen- 
tially independent of frequency implies that the 
waves are deflected at a fixed transverse distance 
from the pulsar line-of-sight by either refraction or 
scattering. This reinforces the conclusion that spe- 
cific isolated structures in the line of sight are re- 
sponsible. The observed frequency independence 
constrains the angular size of a lens, based on the 
law of plasma refraction, to be a small fraction 
of the angle of deflection which we observe up to 
±20 mas. Similarly the size of a scattering center 
would have to be much smaller than this angle of 
deflection. 

The very high apparent axial ratio we observe 
strongly suggests that an ordered magnetic field 
determines the geometry of the scattering. The 
lack of any magneto-ionic signature (i.e., the exact 
match between the dynamic spectra in right and 
left circular polarizations) implies that the radio 
frequency phase imposed by the plasma is simply 
proportional to the column density of the local 
electron density n c . We expect that plasma struc- 
tures will be more tightly confined transverse to a 
magnetic field than parallel to it, i.e., the density 
irregularities will be field-aligned as they are in 
the solar corona. Since we expect such structures 
to scatter preferentially in their narrow dimension, 
our initial idea is that the major axis of the angular 
scattering would be orthogonal to the mean mag- 
netic field. Thus one scenario for the main image is 
a roughly parallel set of filaments (or sheets), con- 
trolled by a magnetic field at right angles to the 
axis of scattering and extending at least over the 
projected length of the image (~ 15 AU). Waves 
passing through a filament are scattered by 10- 
30 mas due to a locally enhanced column density 
of electrons (and presumably higher plasma tur- 
bulence) which would be seen as a single arclet. 
We call this first scenario the orthogonal geome- 
try. The axial ratio R would then correspond to 
the axial ratio of the micro-turbulence (i.e., at the 
diffractive scale). 

There is however, an alternative scenario in 
which a dense filament of plasma is confined by the 
local magnetic field, which is parallel to the axis 
of scattering. We call this the parallel geometry. 



Such a filament must have denser knots of micro- 
turbulence, which are responsible for the discrete 
arclets. The flux ropes seen in the solar wind are 
a possible plasma structure (see Birn et al. (2006) 
and references therein). The denser knots could 
scatter isotropically and so the filament would not 
have to be located exactly in front of the pul- 
sar. In such a case the major axis of the scat- 
tered image (about 3 mas) would be determined 
by diffractive scattering caused by microstructures 
of the order of 3000 km in scale, but the minor 
axis would be set by the thickness of the fila- 
ment. Such scattering has been observ ed in the 
solar corona ( Dennison fc Blessing! 19721 ) and dis- 
cussed in the context of interstellar scattering by 



Cordes fc Laziol ([20011 ) . We assume that the thick- 



ness of such a filament would be of the order of 
the size of a "knot", on the order of 0.05 AU or 
0.1 mas. Thus the filament is at least 16 AU long, 
about 0.05 AU in diameter, but not straight. The 
perpendicular RMS angle of 0.4 mas would corre- 
spond to an RMS irregularity in the filament of 
0.2 AU or 4 times its thickness. Because the scat- 
tering is caused by a very thin structure in this 
model, the electron density in the filament would 
have to be considerably higher than normally ex- 
pected (about 5 cm~ 3 ) to cause the observed an- 
gular scattering. 

In the parallel filament model the 1 ms feature 
is readily explained as a separate offset filament 
of about the same density which is approximately 
parallel to the main filament. However, in the or- 
thogonal geometry the observer would not detect 
the anisotropic scattering from a concentration of 
filaments nearly parallel to those causing the main 
arc. Thus the basic simplicity of this geometry, 
that assumes the scattering is so anisotropic that 
the angles of scattering are transverse to the long 
spatial axis, must be augmented by an indepen- 
dent process of scattering or refraction to explain 
the 1 ms feature. A possible process is a localized 
plasma structure in that part of the sight line that 
refracts the scattered waves towards the observer. 
In this case the angular position of the 1 ms feature 
should have a A 2 wavelength dependence, which is 
inconsistent with that observed as shown in ^3.31 

In either geometry the existence of many sub- 
AU sized discrete structures passing within 15 AU 
of the line of sight to our pulsar raises the more 
general question of their distribution in interstel- 
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lar space. The fact that the structures responsible 
lie at a common distance suggests that the pulsar 
happens to lie behind a single larger region con- 
sisting a cluster of many filaments of plasma or a 
single filament with very comp act conden s ations 

~ (|2005l) 
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(knots). Further the fact that 
observed similar arclets 22 months earlier suggests 
a region bigger than 40 AU. Two ideas that seem 
possible are a shock with multiple nearly parallel 
subshocks seems and the flux rope idea already 
mentioned. The number density of such regions 
can only be crudely constrained by the sparse sta- 
tistical sampling of pulsars with similar multiple 
arclets. We hope to distinguish between these two 
basic geometries by detailed modelling of the prop- 
agation. 

6. Conclusions 

This paper describes a novel VLBI technique 
resulting in a two-dimensional image of the scat- 
tering screen of pulsar B0834+06. The baseband 
data that were recorded allowed high resolution 
dynamic spectra to be produced. The secondary 
spectra produced with the two-hour dynamic spec- 
tra could allow sharply defined arclets to be iden- 
tified with delays as high as 1 ms. 

The scattered image was developed by astro- 
metrically mapping points chosen from the sec- 
ondary spectrum to bright points in the sky plane. 
These points were clustered in two clearly defined 
groups: a primary scattering disk which is elon- 
gated and inclined 27±2° to the pulsar proper mo- 
tion direction and a second, non-colinear, feature 
corresponding to the 1 ms feature of the secondary 
spectrum. Diagnostic measurements place the two 
features at essentially the same distance, 65% of 
the way to the pulsar. The two-dimensional dis- 
tribution of points in the scattered image allows 
both transverse components of the effective veloc- 
ity to be determined via relationships connecting 
the Doppler frequency with location in the image. 

The discrete feature at 1 ms delay contains 
about 4% of the total received power. This fea- 
ture is expected to to be visible only for a few 
weeks during which time its delay should drift as 
the pulsar moves; the impact on timing this pul- 
sar at ~ 327 MHz due to such a feature is a time 
variable wander with magnitude ~ 40 /is. This 
should come as a caution to those aiming to per- 



form precision pulsar timing at low frequencies on 
pulsars that exhibit the extreme forms of scintilla- 
tion that are characteristic of B0834+06. Further, 
pulsars with sub-microsecond structure may ex- 
perience apparent pulse profile evolution yielding 
additional complications in their timing. 

We were able to estimate the effective scin- 
tillation velocity vector, which depends on a 
distance-weigthed sum of the velocities of the 
pulsar, the Earth and the sacttering plasma. 
By using the published proper motion we esti- 
mated the velocity of the scattering plasma to 
be 16 ± 10 km s -1 approximately parallel to the 
scattering axis. Since the errors in this interesting 
result are dominated by the uncertainty in the pul- 
sar proper motion, we have undertaken a new set 
of VLBI measurements to improve its precision. 

The interpretation of ISS in pulsars has of- 
ten assumed isotropy in the scattering. The ex- 
tremely anisotropic scattering found here would 
substantially alter any quantitative modelling of 
the plasma were it to be a common feature in 
other regions on the interstellar medium. A de- 
scription of the underlying plasma physics must 
await a resolution of the two possible geometries 
mentioned in the previous section, but the results 
and the method provide an exciting new glimpse 
of the ionized ISM at scales of 0.1 to 10 AU. 
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